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0\ ■ ABSTRACT 

Q-i| We have modified the standard code for primordial nucleosynthesis to include 

O I the effect of the slight heating of neutrinos by e^ annihilations. There is a 

^ [ small, systematic change in the "^He yield, AY c:^ +1.5 x 10"^, which is 

f^ I insensitive to the value of the baryon-to-photon ratio r] for 10^^° ^ i] '^ 10^^. 

.^ I We also find that the baryon-to-photon ratio decreases by about 0.5% less 

/\ I than the canonical factor of 4/11 because some of the entropy in e^ pairs is 

c^ I transferred to neutrinos. These results are in accord with recent analytical 

estimates. 



1 Introduction 

The concordance between the predictions of primordial nucleosynthesis and 
the observed abundances of D, '^He, ^He, and ^Li is one of the cornerstones 
of the hot big-bang cosmology, and provides its earliest test. Because of this 
and the great interest in the very early history of the Universe, primordial 
nucleosynthesis has been called "the gateway to the early Universe." Further, 
big-bang nucleosynthesis has been exploited to provide the most accurate 
determination of the baryon density |l| and to probe particle physics, e.g., 
the stringent limit to the number of light neutrino species . 

Over the past decade there has been continued scrutiny of primordial 
nucleosynthesis, both on the theoretical side and on the observational side: 
Reaction rates have been updated and the effect of their uncertainties quan- 
tified 0, finite-temperature corrections have been taken into account [^], 
and the effect of inhomogeneities in the baryon density explored [^]; the pri- 
mordial abundance of ^Li has been put on a firm basis 0, the production 
and destruction of D and ^He have been studied carefully 0, and astro- 
physicists now argue about the third significant figure in the primordial ^He 
abundance ^. The result is that the "concordance region" of parameter 
space has continued to shrink. The predicted and measured primordial abun- 
dances agree provided: the baryon-to-photon ratio lies in the narrow interval 
3 X 10"^*^ ^ ''? ^ 5 X 10^^° and the equivalent number of light neutrino species 
Ni, < 3.4 1^. The trend motivates the study of smaller and smaller effects, 
and, in particular, the present examination of the small effect of the heating 
of neutrinos by e^ annihilations. 

To place our work in perspective, let us enumerate the usual assumptions 
underlying primordial nucleosynthesis: (i) Friedmann-Robertson- Walker cos- 
mology; (ii) the input of various nuclear reaction cross sections, the most 
important of which is the matrix element for the processes that interconvert 
neutrons and protons; (iii) N,^ nondegenerate neutrino species, i.e., neutrino 
chemical potentials \^u\ <^ T; and (iv) the complete decoupling of neutrinos 
from the electromagnetic plasma before the entropy in e^ pairs is transferred 
to photons. It is the final assumption that our work addresses. 



It has long been known that neutrino interactions with the electromag- 
netic plasma, e.g., z/ + z/ ^^ e+ + e", i/ + e^ <-> z/ + e^, and so on, become 
ineffective (interaction rate per particle F less than the expansion rate H) 
at a temperature of the order of a few MeV. Since e^ pairs do not dis- 
appear and transfer their entropy to the plasma until a temperature of the 
order of rrie/S ~ 0.1 MeV, one expects that neutrinos do not share in the e^ 
entropy transfer. It then follows that long after the e^ pairs disappear the 
ratio of the photon and neutrino temperatures should be (11/4)^/^ and the 
baryon-to-photon ratio should decrease by 4/11 [p!o[ . 

Neither neutrino decoupling nor the disappearance of e^ pairs are in- 
stantaneous events, and so one might expect neutrinos to share slightly in 
the e^ entropy transfer and to have a higher temperature than (4/ll)^/^T-y. 
Because the yields of primordial nucleosynthesis are very sensitive to the neu- 
tron fraction, which around the time of nucleosynthesis is determined by the 
rate of neutron-proton interconversions through the processes n+e^ ^^ p+i^e, 
n+Ue <-> p+e~, and n <-> p+e~ + z/e, they depend critically upon the neutrino 
temperature as the rates for these processes vary as T^. Even a slight amount 
of neutrino heating is potentially important for the ^He since its abundance 
is now discussed to three significant figures. 

A number of authors have tried to quantify neutrino heating and its ef- 
fect on the ^He yield |Q, |Tl], |T2|, |T3|, |I^ . With the except of the most recent 



work, Refs. |TB|, |T^, previous estimates were "one- zone" calculations, i.e., 
the integrated perturbation to the neutrino energy density 6pi, was calcu- 
lated rather than the perturbations to the neutrino phase-space distribution 
functions. All authors agree that neutrino heating of the electron neutrinos 
increases their energy density by about 1% (slightly less for z/^ and z/^, as 
they only have neutral-current interactions), and with one exception, Ref. 
|l^ , all estimate the change in the mass fraction of "^He synthesized to be of 



the order of Ay ~ 10 ^, though there is no consensus as to the sign of this 



small change; the authors of Ref. |jT2[ estimate the change to be 30 times 
larger, AY ^ -0.003. 

In this work we incorporate the results of the most detailed treatment of 



neutrino heating [Q into the standard big-bang nucleosynthesis code ^5 



For the interesting range of the baryon-to-photon ratio we find a system- 
atic increase in the ^He abundance of AF = 1.5 x 10~^, very close to the 
semi-analytical estimate in Ref. WM; we find similar fractional changes for 



the abundances of the other light elements. By integrating the first law of 
thermodynamics we find that due to neutrino heating the baryon-to-photon 
ratio decreases by about 0.5% less than the canonical factor of 4/11; again, 
in good agreement with the estimate of Ref. [|14 . 



We trace the discrepancy in the predicted sign of the change in ^He yield 
to other authors not considering all of the effects of neutrino heating on the 
"^He yield. In order to check the interesting claim that AY = —0.003 we have 
also modified the nucleosynthesis code to take into account the effect of neu- 



trino heating as computed in Ref. |12]; however, we find that the predicted 



change in He is only AY = +1.1x10 , which is consistent with the more 



detailed treatment of neutrino heating. Since the authors of Ref. [jT2[ give 
few details concerning the changes they made in the nucleosynthesis code, it 
is not possible to explain this discrepancy, though we are very confident that 
the change is not as large as they state. 

Our paper is organized as follows; in the next Section we discuss the 
changes that must be made in the nucleosynthesis code when neutrino heating 
is taken into account and how we implemented them. In the final Section, 
we discuss our numerical results, compare them to previous estimates for the 
change in ^He production, and finish with some concluding remarks. 

2 Modifications to the Standard Code 

2.1 Role of neutrinos 

The slight heating of neutrinos by e^ annihilations causes: (i) small pertur- 
bations to the neutrino phase-space distributions; and (ii) small decrease in 
the temperature of the electromagnetic plasma (at fixed value of the cos- 
mic scale factor R{t)) since neutrinos take energy away from the electro- 
magnetic plasma. To understand how these changes affect the outcome of 
nucleosynthesis, let us first review how neutrinos "participate" in primordial 



nucleosynthesis. 

Neutrinos play several roles; first, in governing the neutron-to-proton 
ratio. Specifically, the electron neutrino and antineutrino phase-space dis- 
tributions affect the rates (per nucleon) for the reactions that interconvert 
neutrons to protons and vice versa, Xnp and Ap„. In the standard treatment 
these rates are computed by integrating the well known tree-level matrix el- 
ement squared over the appropriate (thermal) Fermi-Dirac distributions (see 
e.g., Refs. 0, |10[)- Because of neutrino heating, the electron- neutrino distri- 
bution is given by the usual thermal part plus a small perturbation, which 
results in small changes to the weak rates, 6Xpn and 6\np- 

The other roles neutrinos play involve their contribution to the energy 
density of the Universe. The total energy density determines the expansion 
rate of the Universe: 

... 9 

3 ' ^ ^ 

where ptot = P-y + Pe + Pv + Pb- Because of rapid electromagnetic inter- 
actions the electromagnetic plasma is always in thermal equilibrium so that 
Pem(^7) = Pj + Pe is only a function of the photon temperature Ty. And of 
course, the baryonic contribution to the energy density, ps, is very tiny as 
the Universe at this early time is radiation dominated. 

In the absence of neutrino heating by e^ annihilations the neutrino tem- 
perature just red shifts with the expansion, Tq,^ oc R^^, and the neutrino 
energy density poi, oc l/R^. When neutrino heating is taken into account 

Pu = Pou + hu^ (2) 

where pop is the energy density in all three neutrino species in the absence of 
neutrino heating, and 5py is the sum over all three species of the additional 




energy density due to neutrino heating. In Ref. [|14[ the evolution of the 
perturbation to the phase-space distribution of each species is computed; for 
electron neutrinos 5p/ p approaches about 1.2%, and for /i or r neutrinos 
about 0.6%; thus, bp^j p^-, the average over the three species, is about 0.7%. 
The neutrino energy density also appears in the first law of thermody- 
namics, which governs the rate at which the photon temperature decreases 



with time: 

d[pTOTV] = -pTOTdV, (3) 

where the total pressure ptot = Pem + Pu + Pb, Pem(T^) = P-y + Pe, and 
V = R^. Since we shall assume that the three neutrino species are very 
light, m^ -C 1 MeV, the neutrino pressure p^, = Pu/'i- In the absence of 
neutrino heating the neutrino energy density drops out of Eq. (|^) equation 
since p^y oc R~^. When neutrino heating is taken into account this is no 
longer true; as we shall see, the additional term in this equation involving 5py 
leads to a "back reaction" resulting in a slight cooling of the electromagnetic 
plasma. 

2.2 Alterations 

The integral expressions for the unperturbed weak rates A„p and Ap„ cannot 



be calculated in closed form, but the standard code |]T5| allows for either 
numerical integration at each temperature step, or for the use of a series 
approximation (in 1/T^). We have opted for the numerical routine. The 
perturbations to the weak rates are implemented very simply: the numer- 
ical solutions for 5Xnp and 5\pn calculated in Ref. [l^ are added to the 
unperturbed rates by means of a look-up table. 

The effect of the back reaction of neutrino heating on the electromagnetic 
plasma is more complicated. To begin, it is useful to describe how the evolu- 
tion of the photon temperature is computed. At each time step all the energy 
densities and their derivatives are computed, and then stepped forward in 
time by a Runge-Kutta integrator. The time rate of change of the photon 
temperature can be written as 

dT^ d\n.V dT^ dT^ . . 

Ht^ dt dlnV^ dlnV ^^ 

The first law can be used to calculate dT^/dlnV: 

dT^ Pem + Pem + 4:6p,y/3 



dlnV dpEM/dT^ + d5pu/dT^ 



(5) 



Once the evolution of the photon temperature is known, the evolution of 
all other quantities (light-element abundances and so on) follows as in the 
standard case. For example, the evolution of the baryon-to-photon ratio r) is 
governed by 

d\nr]/dt = -3d\n{RT^)/dt. (6) 

Due to e^ annihilations RT^ is not constant, and r] decreases with time. 

A technical note for the experts; in the nucleosynthesis code the first-law 
expression for dT^/dlnV is actually somewhat more complicated because it 
also takes into account: the slight excess of electrons over positrons (electron 
chemical potential fie of order lO^^^T), the tiny energy density and pres- 
sure associated with baryons, and the bookkeeping associated with nuclear- 
binding energies. Since these effects are small and unaffected by neutrino 
heating, we have left them out of our discussion here. 

3 Results and Conclusions 

The "input data" to the nucleosynthesis code needed to compute the ef- 
fect of neutrino heating on the primordial nucleosynthesis are: 6pu/pu, SXpn, 
and SXnp- We consider two approaches to computing these quantities: (I) 
the detailed Boltzmann treatment where the perturbations to the neutrino 
phase-space distributions are computed [Q; and (II) the bulk heating ap- 



proach, where it is assumed that the distortions to the neutrino distributions 
are thermal and only the bulk transfer of energy from e^ annihilations is 
computed |Tl|, |T2[ . In the bulk-heating approach the effect of neutrino heat- 



ing is a slight increase in the neutrino temperature; we use the results of Ref. 



1^ for ST^^ to compute SXi. While we feel that the first approach is more 
accurate, we have also considered the bulk-heating approach because in Ref. 
[|r^ a very large change in the ^He abundance is claimed, AY = —0.003. We 



refer the reader to Refs. |T^ and |jT^ for details about the two approaches. 



The evolution of the energy transfer from the electromagnetic plasma to 
the neutrinos is shown in Fig. 1 for the two methods of computing neutrino 
heating; asymptotically 6pu/ pou approaches 7 x 10~^. It is heartening that 



these two different treatments agree within 15% or so on the integrated mag- 
nitude of the distortion to the neutrino distributions. One consequence of 
the energy transfer is that there are more electron neutrinos and they have 
higher energies, and so the rates for the processes that interconvert neutrons 
and protons increase. However, there is no free lunch: The temperature 
of the electromagnetic plasma drops since it loses energy to the neutrinos 
[Q. Thus a second consequence of the energy transfer is a decrease in the 



neutron-proton interconversion rates due to the drop in the temperature of 
the electrons and positrons. This is a straightforward — but very important — 
implication of energy conservation. 

The third consequence of the neutrino heating is also related to the drop in 
the temperature of the electromagnetic plasma. At a fixed time, the photon 
temperature is slightly lower than in the absence of heating; equivalently, at a 
fixed photon temperature the Universe is slightly younger than in the absence 
of neutrino heating. As is well appreciated, the ^He abundance is determined 
by the neutron fraction at the onset of nucleosynthesis (Tnuc ~ 0.07 MeV); 
which, in part, is determined by the number of neutrons that have decayed 
by this time. Since the Universe is slightly younger, fewer neutrons will have 
decayed. We dub this the "clock effect." 

To summarize, there are three effects: (i) increase in neutron-proton in- 
terconversion rates due to neutrino heating; (ii) decrease in neutron-proton 
interconversion rates due to the drop in the temperature of the electromag- 
netic plasma; and (iii) the clock effect. The discrepancy over the sign of 
the change in the ^He abundance traces to the fact that with the exception 



of Ref. |Tj], all other authors have only considered the first of these three 
effects. 

First, consider the change due to the distorted electron-neutrino distribu- 
tion. For simplicity let us begin by assuming that the perturbation is thermal 
(method II), characterized by a change in the electron-neutrino temperature 
6T^^ . The change in the neutron fraction X„ at the onset of nucleosynthesis 
due to a change in either the neutrino or the temperature of the electromag- 
netic plasma is found numerically to be \T7\ 



5X^ = -0.15T/T. (7) 



It is easy to understand the sign in Eq. (|^): when the temperature rises, the 
rates for neutron-proton interconversions increase and the neutron fraction 
tracks its equihbrium abundance, X„/(l — X„) = exp(— Am/T), longer, 
which leads to a lower neutron abundance when nucleosynthesis commences. 
What is 5TyJ Since the electron neutrinos have both charged- and 
neutral-current weak interactions, they get more than their share of the en- 
ergy transferred to the neutrinos, about as much as mu and tau neutrinos 
combined. Therefore, 

T,, 4 p., - 8 p, ■ ^ ' 

Fig. 1 shows that ^Pvj Pv ^ 7 x 10^'^, and thus it follows that the change 
in the neutron fraction due to the fact that electron neutrinos are hotter is 
bXl = -2.6 X 10-^ 

This is not the whole story; there is a change in the neutron fraction of 
opposite sign due to the slight decrease in the temperature of electrons and 
positrons, which we also estimate by Eq. (|^). If we assume that electrons 
and positrons are relativistic (a good approximation since the neutron frac- 
tion freezes out at a temperature of about 0.7 Me V) and ignore the small 
differences between e^'s and 7's due to statistics, then electrons, positrons, 
and photons each lose the same amount of energy due to neutrino heating. 
Remembering Sp^u = —^Pu, it follows that 

^ - -i^. (9) 

Comparing Eqs. (H) and (P) we see that the fractional change in the electron 
temperature is —2/3 that of the electron-neutrino temperature, leading to an 
increase in the neutron fraction that is only 2/3 as large, SX^ — 1.7 x 10~^. 
The predicted net change in the neutron fraction is thus 



6X^ = 5x;; + 5X^ ~ -0.1 ( - - - -^^-gx w-\ (lo) 

4/ Pi, 



Jv 



Figure 2 shows (5X„ as a function of temperature; the numerical results agree 
well with this simple analytical prediction. 



Figure 2 also shows the result of incorporating neutrino heating into the 
code via method I, where the distortion is not assumed to be thermal. In 
fact, as discussed in Ref. WM, the perturbation to the neutrino spectra 



is highly nonthermal due to the fact that more high-energy neutrinos are 
produced in the process of neutrino heating since neutrino cross sections 
rise with energy. This excess of high-energy neutrinos further enhances the 
neutron-production rate, which as the temperature drops is becoming more 
suppressed by the neutron-proton mass difference, and therefore we expect 
6X!^ to be larger — which is precisely what is seen in Fig. 2. Since SX!^ is 
larger, the near cancellation between SX!^ and SX^ is even more precise: For 
method I, dX^ ~ -2 x 10"^ 

For reference, the the mass fraction of ^He synthesized is related to the 



neutron fraction at freeze out by: Y ~ 1.33X„ (see e.g., Ref |]T4|)- Thus, 
the predicted change in the ^He mass fraction due to the first two effects is: 
AFi+2 ^ -3 X 10"^ (method I) and -1.1 x 10"^ (method II). 

The clock effect involves the age of the Universe at the epoch at which 
nucleosynthesis commences, T = Tnuc — 0.07 MeV. Since the Universe is 
slightly younger when nucleosynthesis commences when neutrino heating is 
taken into account, fewer neutrons decay from the time that the neutron 
fraction freezes out, leading to a larger ^He abundance. In Ref. [T^ the 
change in the "^He abundance due to the clock effect was estimated to be 
^^ciock — +1-5 X 10~^. Figure 3 shows the total change in ^He abundance 
as computed by our modified version of the standard code. For method I, 
Ay is about +1.5 x 10~^, while for method II it is about +1.1 x 10~^, which 
indicates that the AFciock ~ 2 x 10~^, in reasonable accord with the previous 
estimate. 

There are a couple of fine points to be made about the baryon-to-photon 
ratio. In the standard scenario the baryon-to-photon ratio decreases by a 
factor of 4/11 from its pre- nucleosynthesis value to its post-nucleosynthesis 
value, due to the entropy transfer from e^ pairs to the photons (see Fig. 4). 
When neutrino heating is taken into account the decrease is less, by about 
0.5% (see Fig. 4), which means that for a fixed value of rj today, the value 
of 1] before e^ annihilations was smaller. We remind the reader that one al- 



ways specifies the yields of primordial nucleosynthesis in terms of the present 
value of the baryon-to-photon ratio. This suggests a fourth effect of neu- 
trino heating on nucleosynthesis, involving the fact that the value of rj at 
early times is always smaller when neutrino heating is taken into account; 
we dub this the ?7-effect [|14|. Thankfully, this effect for most values of r] is 



small because somewhat before the onset of nucleosynthesis rj has reached 
its asymptotic (present) value; see Fig. 4. For the most interesting values, 
10^^° < r] < 10~^, Ay is insensitive to the value of r]. For extreme values of 
T] it becomes r] dependent. 

For large values of r], AY depends upon rj because of the ?7-effect: As 
one increases t] the onset of nucleosynthesis occurs earlier; for large enough 
1] it occurs before rj has reached its asymptotic value; thus, the value of 
7] during nucleosynthesis is slightly smaller when neutrino heating is taken 
into account; since Y increases monotonically with t], the amount of ^He 
synthesized decreases due to this effect. This is precisely the behaviour seen 
in Fig. 3: For t] ^ 10^^, AY decreases. 

To understand why AY also decreases for very small values of rj, we must 
first recall why the primordial helium abundance drops so precipitously for 
small values of rj (for t] < 10~^^ the mass fraction of D synthesized is actually 
greater than that of ^He). Small t] means that number densities of all nuclear 
species are small, so that nuclear-reaction rates are correspondingly lower: 
Tnuciear oc Tj. FoT cxtrcmcly low valucs of 1], by the time nucleosynthesis 
commences nuclear reaction rates have been come ineffective (Fnucicar ~ H), 
and the amount of "^He produced depends upon the relative effectiveness 
of the nuclear reactions: Y oc FnuciearZ-f^- Neutrino heating increases the 
expansion rate (at fixed photon temperature), therefore the FnuciearZ-f^ is 
smaller and less ^He is synthesized. This is precisely what is seen in Fig. 3: 
For 7] -C 10~^°, Ay decreases with decreasing t]. 

To conclude, neutrino heating affects the synthesis of ^He in four distinct 
ways; by incorporating the effect of the slight heating of neutrinos by e^ 
annihilations into the standard nucleosynthesis code we have quantified its 
effect on nucleosynthesis and clarified previous conflicting estimates. The 
net result of the four effects is a slight increase in the mass fraction of ^He 
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synthesized, AY ^ +1.5 x 10 , for the interesting range of r]. 
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In Ref. []T4| the back reaction of neutrino heating on the electromag- 



netic plasma is carefully discussed; at fixed value of the cosmic scale 
factor, T-y is lower when neutrino heating is taken into account by 
6T^ = -6p„/{dpEM/dT^) ^ -\{5p^/poi,)T^, which follows by energy 
conservation. 

[17] The empirical relation SXn — —O.lSTi/Ti can be derived. To do 
so, define all perturbations with respect to fixed value of the scale 
factor, cf. Ref. [14|. The freeze-out value of X„ is set by its equi- 



librium value when T/H ~ 1: X„/(l — X„) = exp{—Am/Tf), 
where {T/H)\tj = 1. Write F = aR-^{rl + r^) and H = bR-"^, 
where R~^ = Tq^ (neutrino temperature in absence of e^ heating), 
r^ = T^/Tqj, ~ 1, Tg = T^/Tqi, ~ 1. The effective temperature "felt" 
by nucleons, which determines the equilibrium neutron-to-proton ra- 
tio, Tcff ~ {T,y + T^)/2 = R~^{x,^ + Xe). Since dp-^M = —^Pu neutrino 
heating does not affect the expansion rate. After some algebra it fol- 
lows that 5Xn = -|X„(1 - Xn)iAm/Tf)ST^/Ti ^ -O.lSTi/Ti. 
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FIGURE CAPTIONS 

Figure 1: The evolution of tlie total perturbation to the neutrino density 
due to heating by e^ annihilations as a function of the photon temperature. 



The solid curves are the results from Ref. |IJ]; the broken curves are those 
fromRef. pl. 



Figure 2: (a) The evolution of the neutron fraction X„ as a function of the 
photon temperature: Xn tracks its equilibrium value until T^ ~ 0.3 MeV, 
when it levels off because of the freeze out of the weak interactions; it then 
slowly decreases due to neutron decays; the precipitous drop occurs because 
of the onset of nucleosynthesis (Ty ~ 0.07 MeV). (b) The change in the 
neutron fraction 6Xn due to the effects of neutrino heating as a function of 
the photon temperature: 6Xn begins to level off at Ty ~ 0.2 MeV due to the 
freeze out of the weak interactions; it then rises because at a given value of 
T^ the Universe is younger and fewer neutrons have decayed ("clock effect"); 
it drops to zero when nucleosynthesis commences. The solid curves are based 



upon the results of Ref. [Q; the broken curves upon those of Ref. |T2|; these 
results are for rj = 3 x 10"^°. 

Figure 3: The change in the predicted ^He abundance due to neutrino 
heating as a function of the present baryon-to-photon ratio. The solid curves 



are based upon the results of Ref. IQ; the broken curves upon those of Ref. 
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Figure 4: (a) The evolution of the baryon-to-photon ratio as a function of 
the photon temperature; (b) The in change baryon-to-photon ratio, Ar//?], 
due to neutrino heating as a function of photon temperature. Note, we have 
chosen the initial value of t] with and without neutrino heating so that the 
final value is identical. The solid curves are based upon the results of Ref. 



14]; the broken curves upon those of Ref. 
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